PRODUCTS OF NON-STATIONARY RANDOM MATRICES 
AND MULTIPERIODIC EQUATIONS OF SEVERAL 
SCALING FACTORS 

AI-HUA FAN, BENOiT SAUSSOL, AND JORG SCHMELING 

Abstract. Let f3 > 1 be a real number and M : R — > GL(C d ) be a 
uniformly almost periodic matrix- valued function. We study the asymp- 
totic behavior of the product 

P n (x) = M{j3 n ~ l x) • ■ ■ M((3x)M(x). 

Under some condition we prove a theorem of Furstenberg-Kesten type 
for such products of non-stationary random matrices. Theorems of King- 
man and Oseledec type are also proved. The obtained results are applied 
to multiplicative functions defined by commensurable scaling factors. 
We get a positive answer to a Strichartz conjecture on the asymptotic 
behavior of such multiperiodic functions. The case where is a Pisot- 
Vijayaraghavan number is well studied. 



1. Introduction 
Kingman's subadditive ergodic theorem was originally proved in 1968 



| ECU , Ki2|. A more recent proof was given by Katznelson and Weiss in 1982 
| KW ] . It is one of the most important results in ergodic theory. In this paper 
we consider the following set-up which resembles a dynamical system with- 
out invariant measure and try to get results similar to Kingman's theorem. 
Let (3 > 1 be a positive real number. Let {f n } be a sequence of uniformly 
almost periodic functions (i.e. in the sense of Bohr, see Section |2.1|) defined 



on the real line R. Suppose the following subadditivity condition is fulfilled 

f n+m (x) < f n (x) + f m (f3 n x) for a.e. x and all n,m. 

where a.e. refers to the Lebesgue measure. We would like to study the 
almost everywhere convergence of n _1 / n (x). The Kingman theorem applies 
in the special case where (3 > 1 is an integer and the / n 's are periodic. The 
typical case in our mind is 

/„(*) = log \\M{j3 n - l x) ■ ■ ■ M((3x)M(x)\\ (1.1) 
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where M :R —> GL(C d ) is a matrix-valued uniformly almost periodic func- 
tion. We will prove that the limit lim n _ +00 n~ l f n (x) exists almost everywhere 
(a.e. for short) with respect to the Lebesgue measure under the condition 
that the n~ l f n (x) have joint periods (see Theorem ^^). As a consequence, 
an Oseledec type theorem is proved for the matrix products involved in fll.lD 



(see Theorem 2.£). It is proved that the condition on the existence of joint 



periods is satisfied when (3 is a PV- number (see Section^). 

Our consideration is partially motivated by the study of multiperiodic 
functions, already investigated by Strichartz et al. [ JRS| , Fan and Lau 



and Fan |Fj]. By a Multiperiodic function of one real variable we mean any 
function F : R — ► R which is a solution of a functional equation of the 
following form 

m = fl (±) F (±) + ... + fd (L) F (± 



Pi/ \PlJ \PdJ \Pd 

where d > 1 is an integer; p\ > 1, . . . ,pd > 1 are d real numbers, called 
scaling factors; fi, ■ ■ ■ , fd are d complex valued functions defined on the real 
line, called determining functions. The equation will be called a multiperi- 
odic equation. 

We will assume that the determining functions /,• are periodic or almost 
periodic in the sense of Bohr, as is the case in most applications. We will also 
assume that the scaling factors pj are commensurable in the sense that pj 
are powers of some real number (3 > 1. Without lost of generality, we assume 
that pj = {3 3 for 1 < j < d. Then the multiperiodic equation becomes 

As far as we know, there is few work done for the non-commensurable case 
which is much more difficult. 

In the literature, the case where d = 1 and (3 = 2 (or an arbitrary 
integer) has been studied, especially in the theory of wavelets Q. In fact, 
the scaling function ip of a wavelet satisfies a scaling equation 



<p(x) = y^ j a n ip(2x 



n 



The Fourier transform of <p satisfies a multiperiodic equation of the form 
( |1.2j ) with only one scaling factor (3 = 2 and only one determining function 

h{x) = f{x) = \Y,n a ne inX - 

The scaling functions in wavelets constitute a class of functions sharing 
a kind of similarity. More generally, multiperiodic functions arise as Fourier 
transforms of self-similar objects such as Bernoulli convolution measures 
(d = 1, (3 > 1 being a real number and / being a trigonometric polynomial), 
inhomogeneous Cantor measures (d may be greater than 1) or more general 
self-similar measures produced by iterated function systems (see [g). 

In the case of one scaling factor (i.e. d = 1, then we write f\{x) = f(x)), 



the existence of the solution of the multiperiodic equation (1.2) is simple and 
is assured by the consistency condition /(0) = 1 and a regularity condition, 
say / is Lipschitz continuous. Actually the solution can be written as an 
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infinite product 



00 / \ 



fd 



For the existence of the general equation (1.2), we have 

Theorem A. Let d > 1. Suppose that the determining functions fi, ■ 
are Lipschitz continuous, and satisfy the consistency condition 

/i(0) + --- + / rf (0) = l. 

Suppose furthermore that /i(0), . . . , /d(0) £ [0, +00). T/ien equation (|1.2| 
admits a unique continuous solution F such that F(0) = 1. 

The proof of this theorem is postponed to Section |I~2| 

Our study of equation (fk^) is converted to that of vector valued equa- 
tions of the form 



G(x) 



M 







G 







(1.3) 



where M : R — > A / irf X d(C) is a matrix valued determining function and 
G : R — > R rf is a vector valued unknown function. Matrix products will 
be involved in the study of equation (|l.3[) , which produces some difficulties. 
However, equation (|l,3j ) is a simple recursive relation because it contains 
only one scaling factor. Equation (L2) is equivalent to equation ( |1.3| ) with 
M(x) and G(x) equal respectively to 

/ h(x) / 2 (§) ••• fd-xbfc) fd(^r) \ 



\ 



1 





















1 












/ F{x) \ 



and 



(1.4) 



/ 



We would like to know the asymptotic behavior at infinity of the so- 
lution G. This is a natural question because G often represents Fourier 
transform of a function (a measure or a distribution) and the asymptotic 
behavior at infinity describes quantitatively the regularity of the solution. 
Unfortunately, there is no closed form formula for G in general and the be- 
havior of G is rather complicated, as is shown by the Fourier transform of 
the Cantor measure (d = 1, /3 = 3 and /(£) = cos£). 

Following I JRSfl , we will study the pointwise asymptotic behaviors of 

h n (x) := -\og\F(p n x)\ 
n 

as n — > 00. 

We will prove that, under some conditions, the limit lim n ^oo h n (x) ex- 
ists and is equal to a constant almost everywhere with respect to Lebesgue 
measure. 

This answers partially a questions in || JRSfl . More precisely, we have the 
following results, whose proofs are postponed to Section fO| . 

Theorem B. Suppose that the conditions in Theorem A are satisfied. Fur- 
thermore, suppose that the determining functions /1, • • • , fd o,re either iden- 
tically zero or strictly positive and 1-periodic, and that (3 > 1 is a Pisot 
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number. Let F be the solution of the equation ( |1.2| ) . Then there is a con- 
stant C such that 

lim -F{f3 n x) = £ a.e. 

n— >oo n 

The constant C in the theorem is the leading Liapunov exponent of the 
matrix M(x) above, defined by 

C(M) = inf -M\og\\M(f3 n - 1 x)M(f3 n - 2 x) ■ ■ ■ M(Bx)M(x)\\. 
n>l n 

where M/ denotes the Bohr mean of an almost periodic function / (see 
Section [24] below). 

Theorem C. Suppose that the conditions in Theorem A are satisfied. Fur- 
thermore, suppose that the determining functions /i, • • • , fd o,re 1-periodic, 
Lipschitz, and that (3 > 1 is a Pisot number with maximal conjugate of 
modulus p. Let F be the solution of the equation ( |1.2| ), If 

(1 + \h(x)\ + ■■■ + \h- i (xli3 d - 2 )\){\li(x)\ + -■■ + \fd*/0 d - v )\) , -i 

? < " 

then there is a constant A EM such that 

1 d-l 

lim - log V \F(p n ~ j x) | = A a.e. 

n—>oo n — ' 
3=0 

In this case we do not know if the constant A equals the leading Lyapunov 
exponent of the matrix. 

Theorem A will be proved in Section 4.2 as a special case of a more 
general result (Theorem [4.1|) , Theorm B and Theorem C in Section 4.3. 
Both Theorm B and Theorem C are consequences of our Kingman's Theorem 



(Theorem 2.5 ) and Oseledec's Theorem (Theorem |2.9| ) which are discussed 
in Section 2. In Section 3, we prove that the joint period condition required 
in both Kingman's Theorem and Oseledec Theorem is satisfied when /3 > 1 
is a Pisot number. 



2. Kingman theorem and Oseledec theorem 

2.1. Total Bohr ergodicity and joint e-period. 

Let us first recall the definition of uniformly almost periodic functions and 
some of their properties (see |Bo|]). Next we will introduce the notions of 
total Bohr ergodicity and a joint e-period. 

Let / be a real or complex valued function defined on the real line. 
A number r is called a translation number of / belonging to e > (or an 
e-period) if 

sup|/(x + r) - f(x)\ < e. 

rrGM 

We say that / is a uniformly almost periodic (u.a.p.) function if it is con- 
tinuous and if for any e > the set of its translation numbers belonging to 
e is relatively dense (i.e. there exists a number t > such that any interval 
of length I contains at least one such translation number). H. Bohr proved 
that the space of all u.a.p. functions is a closed sub-algebra of the Banach 
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algebra C&(R) of bounded continuous functions equipped with the uniform 
norm and that it is the closure of the space of all (generalized) trigonometric 
polynomials of the form 

finite 

For any u.a.p. function /, as is proved by Bohr, the following limit exists 



Mf = lim - / f(x)dx. 



It is called the Bohr mean of /. For any locally integrable not necessarily 
u.a.p. /, we define M/ as the limsup instead of the limit. 

Definition 2.1. A sequence of real numbers (u n )n>o is sa id to be totally 
Bohr ergodic if for any arithmetic subsequence (u am+ b) m >o (a > 1,6 > 
being fixed) and for any real p > 0, the sequence (u am+ bx) m >o is uniformly 
distributed (modulo p) for almost every x € K with respect to the Lebesgue 
measure. 

The following is the main property of totally Bohr ergodic sequences that 
we will use. 

Lemma 2.2. Suppose that (u n )n>o is a totally Bohr ergodic sequence. Then 
for any u.a.p. function f and any integers a > 1, 6 > 0, we have 

1 N 

lim — } f(u an+b x) = Mf a.e. 

N^oo iv — 
n=0 

Proof. It is a consequence of the fact that u an +f,x is uniformly distributed ( 
mod p) for almost all x, for any real p > 0, the fact that / can be uniformly 
approximated by trigonometric polynomials and the Weyl criterion. □ 

Remark 2.3. Suppose the (u n ) n >o is a sequence such that inf n ^ m \u n — u m \ > 
0, then the sequence (u n x) is uniformly distributed for almost every point 



x |KLj]. Consequently, the sequence (u n ) is totally Bohr ergodic. A more 
special case is u n = (3 n with f3 > 1. This is the most interesting case for us. 
On the other hand, no bounded sequence can be totally Bohr ergodic. 

Definition 2.4. Let (-F n )n>o be a sequence of u.a.p. functions. Let e > 
and N € N. A real number r is called a joint e-translation number for 

(Fn)n>N if 

sup sup \F n (x + r) — F n (x)\ < e. 

n>N x£R 

If for any e > there exists N(e) such that such joint e-translation numbers 
for (i ? n)n>AT( e ) are relatively dense, we say that (F n ) n >o has joint periods. 

2.2. Kingman's theorem. 

Following ideas of Katznelson and Weiss [KW] we prove the following 



version of Kingman's Theorem. The difficulty in our case is that we have to 
deal with an infinite measure space. We are also dealing with non-stationary 
sequences. 
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Theorem 2.5. Let (it n ) n >o be a totally Bohr ergodic sequence of real num- 
bers and (/ n )n>o be a sequence of uniformly almost periodic functions. Sup- 
pose 

(i) The sequence (n~ 1 f n ) has joint periods. 

(ii) The following subadditivity is fulfilled 

fn+m(x) < f n (x) + f m (u n x) for a.e. x and all n,m. 

(iii) For any n > 1 

sup (f m (u n x) - f m {x)) < oo for a.e. x. (2.1) 

m 

Then the following limit exists and is a constant 

1 1 
lim —f n (x) = inf — M/ n for a.e. x. 

n-^oo n n n 



Proof. The proof is a modification of Katznelson- Weiss' proof |KW], With- 
out loss of generality we assume that uq = 1. 
Let 7 = inf n -M.f n . Let us put 

f~(x) = lim inf -f n (x), f + (x) = lim sup -f n (x). 

n-*oo n n->oo U 

We remark that the subadditivity implies that f (x) < f^{u n x) for all 
n € N and a.e. x € R. In a finite measure space this would imply the 
invariance a.e. In our case it is the boundedness (^^) which implies that 
f^{u n x) < / ± (x) a.e., what makes the functions f ± invariant in the sense 
that / ± (x) = / ± (ti ra x) a.e. (Vn). 

The first part of the proof, i.e. / + < 7 a.e., is simple. We just exploit the 
fact that any (infinite) arithmetical subsequence of (u n x) is Bohr-uniform 
distributed. This provides us with a kind of ergodic theorem. In fact, fix an 
integer N. For any integer n write n = mN + r with < r < N. We have 

m—l 

fn(x) < fN(u k Nx) + fr(u mN x). 

k=0 



The functions / r (r = 0, 1, • ■ • , N — 1) being bounded, by Lemma 2.2 
this readily implies by the Bohr- uniform distribution of (ukNx)k>o that 

1 m_1 1 
f+(x) < lim — V f N (u kN (x)) = -Mf N 
m^oo rnJS + r z — ' Jy 

for a.e. x. Hence f + (x) < 7 for a.e. x. 

Next, we want to prove that f~(x) > 7 for a.e. x. For this we assume 
that 7 > —00, otherwise it is trivially true. Adding to each f n the constant 
value — n||/i|| creates a new subadditive sequence f n := f n — n\\fi\\ with 
fn < (n > 1), /- = /- + H/ill and 7/ = 7^ + ||/i||. So, we may assume 
that f n < (n > 1). We can furthermore set /1 = (observe that this will 
not affect the subbadditivity condition since f n < 0). Then for any A > 
we truncate the function /„ in the way 

fn.A = max(/ n , -nA). 
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Note that the sequence f n ,A fulfills the assumptions of the theorem. Note 
that in this case 7 a > —A and also f^(x) > —A for all x. It is clear that 

/±(x)=max(/ ± ( 3 ;),-A). 

Assume that we proved the theorem for the sequence f n ,A for any A. Then 
we claim that /a 0*0 \ f~(x) for a.e. x as A goes to oo. In fact, if 
/a 0*0 > — A for some A then f^(x) = f~(x) and 7a = 7/ > — oo. On 
the other hand if /a 0*0 = — A for all A then obviously f~{x) = — oo and 
7A = 7/ = — oo. This proves the theorem for the sequence f n . 

^From now on we assume that f\ = 0, f n < and f n is truncated and 
we skip the subscript A. Let e > 0. By the hypothesis (i) on the joint 
periods, there is an integer N(e) such that the joint e-translation numbers 
are relatively dense. For these numbers r we have 

f n (x + r) f n (x) 



n n 



<e (\/x S R,Vn > N(e)) (2.2) 



Notice that there is no loss of generality to suppose that N(e) increases as 
e decreases to 0. We define 

nAx) = min \ n > N(e) : —f n (x) < f~(x) + e 
[ n 

Let A e K = {x: n e (x) > K}. Notice that if e' < e" , we have N(e') > N(e") 
and rv 0*0 - n e"{x), so that A e £ C A e % for K > JV(e'). 

We claim that 
Lemma 2.6. For any e > 0, we have 

lim M(A e K ) = 

where WIA = Ml^ denotes the Bohr mean of the characteristic function of 
the set A (defined if necessary with the limsup). 

In order to prove this Lemma ^1] we need the following lemma which says 
that A e K is to some extent periodic. 

Lemma 2.7. For any joint e-translation number r of (/n A**)n>iV(e) > we 
have A 2 £ + t C A e K . 

Let us first prove Lemma |2.7| Suppose x G A 2 £ + r, i.e. x — r € Aj^, then 

fn(x-r) > r ^_ r - ) + 2e (N(2e) <n<K). 
n 

This, together with the fact that r is a joint e-translation number for all 
f n /n with n > N(e)(> N(2e)) (see ([DP), implies 

Ux) > fn(x - r) _ £ > _ + e (7V(e)<n<K). 
n n 

That means x G Thus we have finished the proof of Lemma |2.7| . 

Now let us prove Lemma |2.6| . Since joint ^-translation numbers are rel- 
atively dense there exists L = L(|) > such that any interval of length L 

contains such a joint |-translation number. Since HkA^ = 0, for any r\ > 
there exists Kq > such that 

Alf][-L,L] <L V (yK>K ) (2.3) 
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where | • | denotes the Lebesgue measure (see the definition of n £ (x)). We 
claim that M(A e K ) < r/ (V-RT > K ). Otherwise M(A e K ) > r] for some 
K > Kq. Then by the definition of M.(A e K ) there exists xq G R such that 



xo+L 



XA'(x)dx > Li]. 



xo 



Take a joint ^-translation number r S [— xq — L, —xo], i.e. — L < xq + t < 0. 
Then by Lemma [2.7|, we have 



Alf][-L,L] 



> 



[•XQ + T- 


YL 

7 


Ixq+t 




rxo+L 






X A 


rxo+L 


XA 







X.&(x)dx 



> 



For the first inequality we have used the fact that [xo+t, xq+t+L] C [-L, L] 
and for the last inequality we have used Lemma 2.7. What we have deduced 
contradicts (|2.3|), Thus Lemma |2.6| is proved. 

We continue our proof of Theorem |2.5| . Let S := ||/ _ 
such that M(A e K ) < e/S. We define first 



< oo . Let K be 



f-{x) ifx#A e K , , . 
g{x) = ^ n _ r A f and m{x) 



n e (x) ifx£A e K 







if x G A 



A' 



if x G A 



K 



(2.4) 



(2.5) 



Lemma |2.6| implies that 

Mg <Mf~ + e a.e.. 
Moreover by the invariance of /~ we have (remember that uq = 1) 

g(x) < g(ukx) for a.e. x and all < k < m(x) — 1. 
Then we have 

m(x)—l 

fm(x)(x) < (g{x) + e)m(x) < 9{ukx) + em(x). 

We define inductively mo(x) = and 

m k (x) = m k -i(x) + m(u mk _ i{x) x). 



Now choose R > K and let k(x) be the maximal k for which rrik{x) < R. 
Note that R—m^r x \{x) < K. Now we get by the subadditivity and Equation 



(2.6) 
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fn(x)< fm(u mk{x) x){ u m k {x) x ) + f R-m k(x) (x)( u m k(x) (x) x ) 

k=0 " « ' 

<0 

k(x)—l m k (x) 

< ^2 ^2 9{ujx) + (m k (x) - m k „i(x))e 

k=0 j=m k _ 1 (x) 
m k{x ){x) 

^ ^2 9(ujx) + m k{x) (x)e 
3=0 
R-l 

< 9(ujx) +Re + KS. 

Taking the Bohr mean, using that Mg = M(g o uj) and dividing by R gives 

1 KS 
-Mf R <Mg + e+ — 

KS 

<M/- + 2e+ — 

by Equation ( |2,4[) . Now we let R — > oo and we get 

7 < M/~. 

We claim that /~" < 7 implies /~ = 7 for a.e. x. Suppose this was not the 
case, then one could find e > 0, 5 > and an interval J = (0, L) of length 
I J| = L = L t such that \An J\ = 5 > 0, where 

A = {x: f~(x) <7-e} 
By the invariance of /~ we have u k A = A for all k £ N. Hence, for all k £ N 

1 r feL 1 /" L 5 

XA dx = — / xa ctx > — > 0. 



^fc-^ jo ^ Jo L 

Since limsup fc Ufc = +00 (see remark |2~lj| ) , we have Mi > and thus 

Mr <7(l-v) +(7-e)| <7- 



□ 



Remark 2.8. One can prove a similar theorem for more general sequences 
In this case it seems to be necessary to assume L 1 Bohr-uniform 
distribution. 

2.3. Oseledec theorem. 

Kingman's theorem implies the following Oseledec type theorem (see Ru- 
elle 



Theorem 2.9. Let j3 > 1 be a real number. Let M: R — > GL^C) be a 
uniformly almost periodic function. Write 

M™ = M{f3 n ~ 1 x) ■ ■ ■ M(f3x)M(x). 
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Suppose the q-exterior products ^ log ||(M™) A9 || have joint periods, for q = 
1, . . . , d. Then there is V C M with j3T C T of full Lebesgue measure (in the 
sense that M\T has measure) such that if x £ T then 

a) lim^oo (M™M%) ~t = A x exists. 

b) Let expA^ < • • • < exp A^ be the eigenvalues of A x (where s = s(x) 
and the X% are reals), and . . . , TlP the corresponding eigenspaces. Let 
rrix = dim llP . We have XlJ. = X% and rrP = rrp and 

lim - log ||M>|| = A^ when v £ VP \ V^ r ~ l) 

n—>oo n 

forr = 0,...,s where vP = {0} and vP = U [ x 1] + • • • + UP . 

c) Moreover vP depends measurably on x and M x vP = V^J. 

(r) (r) 

d) In addition the functions X x and m x are constant a.e. 

Proof. This theorem follows in the standard way from the a.e. convergence 
of 

lim -log\\(M^) Aq \\ 

n— >oo n 

for 1 < q < d, which in the classical case is insured by Kingman's theorem. 
So we only need to check that for the functions fP(x) = log || (M™) Ag || and 
the sequence u n = (3 n the assumptions of Theorem 2J5 hold. 

First we note that M _1 is uniformly almost periodic because M is uni- 
formly almost periodic and M^ 1 (x) € GL^(C). Second we note that M Aq 
and (M hq )~ l are again uniformly almost periodic, since each entry is a 
rational function of the entries of M and M~ 1 , respectively. Hence, 

sup ||M(x) A l + sup \\(M Aq )- 1 (x)\\ = W q < DO. 

Subadditivity is obviously fulfilled since (M£+ m ) A 9 = (M^ x ) Aq (M^) Aq . 
Condition (2A) follows from 

r» \^\\ = ||M(/3 n x) A9 (M^) A9 (M(x) A9 )- 1 || < ||(M") A9 || + W q . 

Finally by Remark the sequence f3 n is totally Bohr ergodic, so Theo- 
rem 2J3 applies. Assertions (a), (b) and (c) follows from Proposition 1.3 



(see also the proof of Theorem 1.6) in fRuf . 

Now we prove d). By Kingman's theorem (Theorem |2.5| ), we have for 
almost all x 

lim -log||M" Afc || = inf -Mlog ||M" Afc || (1 < k < s). (2.1) 

n^oo n n>l n 

On the other hand, by the properties of exterior product, if we write k = 
Xyl=i m x + ^ with < j < s and < £ < m x s we have 

lim - log ||M; Afc || = m x s hi s) + ■■■ + mi s - j+1 h ( x s - j+1) + £X x s ~ j \ (2.2) 

n— >oo n 

We can solve Ai" (1 < r < s) from the system ( |2 . 1| ) - ( pT2[ ) . The solution is 
independent of x since it depends only on the right hand side terms in fl2.1|). 



(r) 

Consequently m x is also independent of x. □ 
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3. When /? is a Pisot-Vijayaraghavan number 

We restrict our attention to the special case where j3 > 1 is a Pisot- 
Vijayaraghavan (PV) number and f n are defined by We will prove 

that, under some extra condition, the sequence n _1 / n has joint periods and 
the Kingman theorem and the Oseledec theorem apply. To do this, we need 
a distortion lemma and some properties of PV-numbers. 

3.1. Distortion lemmas. 

Lemma 3.1. Let M: R — > GL d (C) such that 

D := sup ||Af(x)|| ||M(ar) _1 || < oo (3.3) 

Let (x k ) and (y k ) be two sequences in M and let 9 k = \\M(x k ) — M(yf.)\\. 
Then for all n £ N and ^ v € C d we have 



\M{x 1 )---M{x n - 1 )M{x n )w\ 



<exp \C^2D% . 



k=\ 



\M{y l )---M{y n - l )M{y n )w\ 

where C = (sup^^ ||M(x)||) -1 . 
Proof. Let 

Qi(x) = I, Q k (x) = M(xi)M(x 2 )---M(x k _ 1 ) (1< k < n) 
Q n (x) = I, Q fc (x) = M{x k+l )M(x k+2 )---M{x n ) (l<k<n). 
We can write 

IQn(a)v| = tt lQ fc (?/)M(x fc )Q fc (x)v| 

IQn(y)v| 1 = 1 |Q fc (y)M(y fc )Q fc (x)v| • 1 • J 

Setting E k = M(x k ) ■ M~ 1 (y k ) — I we get the following estimate for the 
numerator of the general term in the above product 

\Q k (y)M(x k )Q k (x)v\ = \Q k {x){l + E k )M{y k )Q k {x)v\ 

= \Q k (x)(I + E k )Q k (x)- l Q k (x)M(y k )Q k {x)v\ 
= \(I + E k )Q k (x)M(y k )Q k (x)v\ 

where E k = Q k {x)E k Q k (x) . It follows that 

\Q k (y)M(x k )Q k (x)v\ ~ 

\Q k (y)M(y k )QH^\- l{I + Ekl (3 - 5) 

\\E k \\ <D k - 1 \\E k \\. (3.6) 



It is obvious that 
On the other hand 



\\E k \\ < supllAf^)" 1 !!^. (3.7) 

X 

By combining (|3^) , JjOD , flO| ) and (p?7D, we obtain 
|M(xi)---M(x„_i)M(x n )v| ""' 



|M(yi) • • • MCi/n-ijAfd/^vl A (j 



□ 
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If M(x) is non-negative, the next lemma shows that condition ( |3.3| ) is not 
needed for positive vectors 

Lemma 3.2. Let M : M —* GL^(]R) 6e suc/t i/ioi i/ie entries of M(x) are 
either identically zero or bounded from below by a positive number 5 > 
(independent of entries). Then for any sequences and (y^) in R and 
for any non-negative vector v we have 

|M(ai)---M(a; ra -i)M(a w )v| < 
|M(?/ 1 )---M( 2 / n _ 1 )M(y n )v| " eXP 

where 9k = \\M(xk) — M(yk)\\ and the norm \v\ = Xa=i 072 ^ * s chosen. 
Proof. We may write 

|M(si)---M(x n _i)M(a; n )v| 
^ Af(xi) io ,i 1 Af(a!2)i 1 ,i a ---M(a? n )< n _ 1 , in Ui„ 

We have a similar expression for |M(yi) • • • M(y n _i)M(y n )v|. Now com- 
pare the two expressions term by term. By the hypothesis, both quantities 
A^(^l)to,ii an d ^(yi)i'o,ii are either zero or larger than 5. So, using the 
trivial inequality x/y < e x ^ y ^ 1 we have 



M(xi) io>il < M(yi)j 0jil e 



The same estimates hold for other pairs M(xk)i k _ 1: i k and M{y] t )i k _ 1 ^ k . The 
desired inequality follows. □ 

3.2. Two properties of PV numbers. 

Let j3 > 1 be a PV-number of order r. We denote its conjugates by 
(3[, ■ ■ ■ , P'r-x- Then for ra > 1, denote 

F n = (3 n + (3' n 1 + 
Lemma 3.3. T/ie number F n is an integer and we have 

\P n -F n \ < (r-l)p n (Vn>l) 

where p = max\<j< r -i < 1. 

Given any real number /3 > 1 (not necessarily integral), we can expand 



each number a; € [0, 1) in a canonical way into its fj-expansion [Re] (see also 
II): 

y~y e n (x) 



n=l 



where (e n (cc)) n >i is a uniquely determined sequence in {0, 1, • • • , [/3]} N . We 
may also call (e n (x)) n >i the /3-expansion of x. We note that not all sequences 
in {0, 1, • • • ,[/3]} N are /3-expansions. Let Dp be the set of all possible f3- 
expansions of numbers in [0, 1). A finite sequence e = (ei, ■ • ■ , e n ) (of length 
n) in {0, 1, • • • , [P]} n is said to be admissible if it is the prefix of the f3- 
expansion of some number x. For such an admissible sequence, we define 

I(ei, • ■ ■ , e n ) = {x € [0, 1) : e x (x) = ei, • • • , e n (x) = e n }. 
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It is known that if Dp is endowed with the lexicographical order, the map 
which associates x to its /3-expansion is strictly increasing. The set /(e) = 
) is an interval, called a j3 -interval of level n. Its length is denoted 

by |/(e)|. 

Lemma 3.4. Suppose (3 > 1 is a PV-number. There is a constant C > 
such that 

C~ l r n < |/(ei,-" ,e n )\ <Cf3~ n 
for any integer n > 1 and any (3-interval I(ei, • • • , e n ). 

See [H for proofs of Lemma 3.3 and Lemma I 



3.3. Existence of joint periods. 

Definition 3.5. Let j3 > 1 be a positive real number and let M : R — ► 
7V4,ixd(C). If the entries of M are functions of the form f(f3 n x) (n € Z) 
where / is 1-periodic continuous, we say that M is (3-adapted u.a.p. 

Remark 3.6. The matrix M(x) defined by ( |1.4j) associated to a multiperiodic 
function is /3-adapted. 

Proposition 3.7. Zei /? > 1 be a PV-number and let M :M.—> GL d (C) be 
(3-adapted and a-Holder continuous. Suppose that 

Dp a < 1 (3.8) 

where p is the maximal modulus of the conjugates of (3 and D is the same as 
in the distortion lemma (Lemma \3.1\ ). Then for any 1 < q < d the sequence 

( \ 

n~ 1 f n (x) has joint periods, where 

fto)( x ) = log \\(M(/3 n ~ l x) ■ ■ ■ M((3x)M(x))) Aq \\. 

Proof. Since M is /3-adapted, the entries of M{j3 k x) are all of the form 
hij(/3 ei 'ix) with 1-periodic function hij(x) and integer i^j > for suffi- 
ciently large k. So, if necessary, we consider log \\{M((3 n ~ 1 x) ■ ■ ■ M(/3 fc °x)) Al3 || 
for some sufficiently large but fixed ko > 0. 

Consider r = (3 m "q m + • • • + (3r]i + rjo where m > 1 and < ry, < (5 
are integers. We are going to show that all such r are joint e-translation 
numbers for n^ 1 f n n \x) with n > N(e), where N(e) depending on e is an 
integer to be determined. 



By Lemma 3.3, we have 



inf \f3 k r-j\ < C'p k 

for all k and some constant C independent of k and r. For k > each 
entry of M Aq (f3 k+k " x) is a degree q polynomial in d 2 variables of the form 
h{(5 i+k x), with h a-H61der, 1-periodic, and t > 0. Notice that we have 
h{/3 e+k (x + t)) = h{(3 l+k x) + 0{p ka ), hence 

\\M{f3 k+ko (x + r)) A9 - M{(3 k+ko {x)) Aq \\ = C qP ka 
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for some constant C q . By the distortion Lemma 3.1 and the above estimate 
we have 



\f, 



(9)(x + r)-/W(x)| 
\\(M(/3 n - 1 x + (3' 



QTl— 1, 



log 



M{(3x + (3t)M{x + t)) a Q\ 



\\(M((3 n - l x) ■ ■ ■ M{[3x)M{x))^i\ 

: C 



c q c'j2 Dk p (k ~ ko)a <^' Dp k ° 

k=l 



So, we may choose N(e) = C/e. In order to finish the proof, it suffices to 
notice that Lemma [3.4| implies that all these r form a subset with bounded 
gap in R. □ 
Proposition 3.8. Let [5 > 1 be a PV~number and let M : R -> GL d (C) be 
f3-adapted and a-Holder continuous. Suppose that the entries of M(x) are 
either identically zero or larger than a constant 5 > 0. Then n f n {x) has 
joint periods, where 

= log \\(M{[3 n - l x) ■ ■ ■ M(Px)M(x))^)\\. 

Proof. The proof is the same as the last proposition. But we use Lemma |3.2| 
instead of Lemma |3.l|. □ 



4. MULTIPERIODIC FUNCTIONS 

As we pointed out in the introduction and as we will see in Section [4.2| , 
our scalar equation ( |1.2| ) can be converted to the vector equation (|1.3|). So, 
we first study the vector equation (11.31). 



4.1. Equation G{x) = M(x/P)G(x/P). 

Let M : R — > Mdxdi^) be a matrix valued function. We consider the 
following vector valued equation 

G(x)=m( |W / ' 



3) \P. 

where the unknown G : R — > C d is a vector valued function. 

Theorem 4.1. Lei /? > 1 6e a real number and M(x) be a complex matrix 
valued function. Suppose that M is Lipschitzian and that M(0) is non neg- 
ative and has 1 as a simple eigenvalue with a corresponding strictly positive 
eigenvector v. Then there exists, up to a multiplicative constant, a unique 
continuous solution G(0) ^ of the equation G(x) = M(x/ (3)G{x/ (3). The 
solution can be defined by 



where the convergence is uniform on every compact subset in R. 
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Proof. Write v = (v±, ■ ■ ■ , Vd) 1 ■ We introduce the following norm for C d 

\\z\\ = max — (z = (zj)i <j<d G C d ). 

l<j<d Vj 

Then a matrix A = (dij) G Mdxd(C), considered as an operator on the 
normed space (C d , \\ ■ ||), admits its operator norm 



1 d 

\A\\ = max — > la, 

Ki<d Vj f 



Notice that ||M(0)|| = 1 because M(0)v = v. 

Since the eigenvalue 1 of M(0) is simple (and isolated), and M(x) 
is Lipschitz continuous, by the perturbation theory of matrices, there is a 
neighborhood of 0, say [—5,(5] (5 > 0), such that for any x G [—5,5], M(x) 
has a simple eigenvalue \{x) and a corresponding eigenvector v{x) satisfying 

\\(x)-l\<C\x\, \\v(x) - v|| < C\x\ (x G [-5,5]) (4.1) 
for some constant C > 0. We claim that the limit 

G(x) = lim M ( M ( -^V • • M ( v (4.2) 



exists (uniformly on any compact set). It is clear that the limit function is 
a solution. 
Denote 

Q »w = M (!) M (lO'" M (f; 

The proof of the existence of the limit in is based on the following 

lemma. 

Lemma 4.2. For any 5 > 0, there exists a constant D > such that for 
any n > 1 and any x G [—5, 5] we /iaue 

||Q„0e)||<£> \\Q n {x)w-v\\<D\x\. 

To get the boundedness of |j<3 n ( 3; ) ||> h suffices to notice that 



HQn(x)||<n/ 
i=i 



where the scalar function /(x) = ||M(x)|| is Lipschitzian and /(0) = 1 (we 
have used our choice of the norm of C d ), and that the products converge 
uniformly on [—5,5] to a continuous function [FL|. Now we prove that 

||Q n (a0v-Q„-i(a0v|| <C'M (4.3) 

where C > is some constant. In fact, since M(x)v{x) = \(x)v(x), we 
have 
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Multiplying both sides by Q n -i(x), we get 

\\Q n (x)v - Q n -i(x)v\\ 
x 



< 

Notice that 



Q n {x) v-d 



+ 



Qn-l(x) V - A 



||A(x)v(x) — v || < \\X(x) — 1|| 11^(2;) || + \\v(x) — v||. 



Using the last inequality, the estimates in fl4.l|) and that we have just proved 
IIQnX^)!! < D, we obtain (f4.3|) . Then for n > m 

-A c'\x\ 

\\Q n (x)v - Q m (x)v\\ < 22 \\Qk(x)v - Q k -i(x)v\\ < 

fc=m+l 



That means Q n (x)v is a Cauchy sequence in the space C([—S,5]) of contin- 
uous functions equipped with uniform norm. Since for any fixed integer no, 
we have 



lim Q n {x)v = Q no {x) ■ lim Q n — v, 

n— >oo n— »oo \ p' l J 

it follows that the uniform convergence of Q n (x)v on [—5,5] implies its 
uniform convergence on any compact set. 

The uniqueness of solution is easy. Let G / be a solution. First notice 
that G(0) is an eigenvector of M(0) associated to the simple eigenvalue 1. 
Hence we may assume that G(0) = v. By iterating the equation, we get 

G(x) = Q n (x)G (jL\ = Q n {x)v + Q n (x) (g - 

The last term converges to zero (uniformly on any compact set) because of 
||Qn(^)|| — D. So, G{x) must be the limit of Q n (x)v. □ 



Remark 4.3. In the theorem, neither the almost periodicity of M{x) nor the 
positivity of M(x) is required, but only the positivity of M(0). That 1 is 
an eigenvalue of M(0) is necessary for the equation (|1.3| ) to have a solution 
G(x) such that G(0) + 0. 

Remark 4.4. The Lipschitz continuity is not really necessary. Holder conti- 
nuity or even Dini continuity is sufficient. 

Remark 4.5. If the entries of M are (real) analytic, then the solution G 
is also analytic. Because, for any xq € R, there is a disk on the complex 
plane centered at xq on which Q n (x)v (as functions of complex variable x) 
uniformly converges. 

4.2. Existence of multiperiodic functions. 

Here we give a proof of Theorem A based on Theorem 4.1. 
Let M(x) be as in (1.4). It is easy to see that the characteristic polynomial 
of M(0) takes the form 

P(u) =u d - h(<S)u d - x - f 2 (0)u d - 2 /d_i(0)« - fJO) 
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The consistency condition implies that 1 is an eigenvalue of M(0). Notice 
that 

P\l) = /x(0) + 2/ 2 (0) ■■■ + (d- l)/ d -i(0) + d/ d (0) > 0. 

So, the eigenvalue 1 is simple. By Theorem |4.1| , there is a unique solution 
of G(x) = M(x/0)G(x/P). Let 

G(x) = (G l (x),G 2 (x),--- ,G d {x)) 1 

Then G\(x) is a solution of ( |1.2| ). If i 7 is a solution of ( |1.2| ). Let 

Gi(x) = F(x), G 2 (x) = F(x//3), • • • , G d {x) = F{x/f3 d - 1 ). 

Then G = (G lt • • • , is a solution of = M(x / (3)G(x / f3). Thus the 
uniqueness of the solution of Equation ( |1.3| ) implies that of Equation ( |1.2|) . 

4.3. Asymptotic behavior of mult iper iodic functions. 

Let us consider the asymptotic behavior of a multiperiodic function, or 
more generally the asymptotic behavior of a solution G of Equation ( |1.3D 
provided it exists (the existence may be guaranteed by Theorem |4.1| . 

Theorem 4.6. Let (3 > 1 be a PV-number whose maximal conjugate has 
modulus p. Let M : M. — > GL^(C) be a (3-adapted u.a.p. Holder function 
of order a > 0. Suppose that G is a solution of G(x) = M(x/(3)G(x/f3). 
Suppose furthermore that one of the following conditions is satisfied 

(i) Dp a < 1 where D = sup xm |] Af (a:) |[ ||iWf (a:) -1 |J (NB. (3 must be Pisot). 

(ii) The entries of M(x) are either identically zero or larger than a con- 
stant 5 > 0. 

Then for a.e. x € M the limit 

h(x) = lim -\og\G((3 n x)\ 

n—foo n 

exists and is independent of x. 



Proof. We first consider the case (i). By Proposition 3/7, Theorem ^9 ap- 
plies. Hence for a.e. x, if we denote by r{x) the integer such that the vector 

G{x) G VP \ yj r_1) we get 

lim -\og\G{(3 n x)\ = lim - log \M%G{x)\ = 

n— >oo n n— >oo n 

But G{f3x) = M{x)G{x) hence G(f3x) G vjg \ V^ x ' 1] , from what follows 
r(j3x) = r(x), i.e. r is invariant. Hence constant a.e. because of the total 
Bohr ergodicity of the sequence P n . 



Case (ii). We use the notation of Proposition Since n 1 f n has joint 



periods, Theorem 2.5 applies (see the proof of Theorem 2.6 for details). 



Hence the following limit exists a.e. 

lim -f n (x) = C, 

n— >oo n 

where C = inf n ±M(/ n ). In view of G{f3 n x) = M(/3 n ~ 1 x) • • • M{x)G{x) the 
positivity of M and G gives 

lim - \og\G{f3 n x)\ = C a.e. 

n— »oo n 
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□ 

Note that when G is the solution of equation (|l.3j ) with M and G given 
by (1.4) we have 

\G{x)\ = \F{x)\ + \F(x/(3)\ + ■■■ + \F(x/(3 d - l )l 



(4.4) 



thus the asymptotic behavior of ^ log \G{f3 n x)\ and -log ^7=0 \F{(3 n ~i x)\ 
are the same. Thus Theorem C follows as an immediate corollary of The- 
orem [0| This partially answers a question in [ JRS| (Conjecture 4.1., p. 
263). 

We prove now Theorem B. By the primitivity of M(0) and the hypothesis, 
there exists an ingeter t > 1 such that M(x) := M(x/(3 T ^ 1 ) ■ ■ ■ M{x/j3)M{x) 
has all its entries strictly positive (even larger than c5 T for some constant 
c > 0). Consider the equation 

G(x) = M(x/f3)G(x/fi T ). 

We examine the first entries of both sides. We can find two constants < 
C\ < C2 such that we get 

Cl F(x) < F{x/(3 T+1 ) + ■■■ + F{x/f3 T+d ) < c 2 F(x). 

Thus Theorem B follows from Theorem 



Theorem 4.7. Under the same conditions as Theorem \4- o\ , for any q € M. + , 
the following limit exists 



lim — log 

n^oo TL 



\\M(f3 n - l x) ■ ■ ■ M{(3x)M{x)\\ q dx. 



Proof. Write 



P n {x) q dx with P n {x) = \\M(p n ~ l x) ■ ■ ■ M{(5x)M{x 



o 



It suffices to show that there is a constant C > such that 



Z, 



n+m 



< CZ n Z m (n>l,m> 1). 



We assume that q = 1, just for simplicity. We will use the fact that there is 
a constant L > such that 



\\M(x) - M(y)\\ 2 < L\x -y\ a (Vx,yGM). 

M n M n _i ■ ■ ■ MiMq for the (noncommuta- 
tive) product of the matrices M , . . . , M n . Write 



We use the notation Tlf =0 M,, 



J n+m 



< 



E 



E 



m— 1 n—1 
H M{(3 n+k x) ■ Y[ M(/3 j 

k=0 j=0 



dx 



n—1 



n M{p 

3=0 



m—1 



H M(l3 n+k x) 



k=0 



dx 



where the sum is taken over all /3-intervals /(e) of level n (see Lemma 3.4). 
Let a e be the left endpoint of /(e). The integral in the last sum , after the 
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change of variables f3 n (x — a e ) = y, becomes 



n \I(e)\ 



n-1 



HM(pia e + f3- n+j y) 

3=0 

Notice that if e = (ei, • • • , e n ), then 



m—1 



[] M{f3 k y + f3 n+k a e 



k=0 



dy. 



(3 n+k a e = P n+k (^ + ■ ■ ■ + = f3 n+k ~ l e 1 + ■■■+ (3 k e n . 



So, by Lemma 3.2, there is an integer n e such that 

\p n+k a e - n e \ = 0(p k + p k+l + ■■■ + p 71 ^- 1 ) 



0(p k ) 



(4.5) 



By the distortion lemma (Lemma 3A, Lemma 3.2), we have 



n-l 




n— 1 


l[M(pa £ + p- n+j y) 


< c 




3=0 




3=0 


m—1 

[] M{(3 k y + (3 n+k a € ) 

k=0 


< c 


m—1 

II M (P k y) 

k=0 



Therefore, we get 



Z n +m < C(3 n ^^P n (a e ) < C'Z n Z n 



□ 



Corollary 4.8. Let F be the multiperiodic function defined by ( |1.2| ). Sup- 
pose that (3 > 1 is a PV-number and that fi,--- , fd are either identi- 
cally zero or larger than a constant 5 > 0. Suppose further that M%, := 
M((3 f ~~ 1 x) ■ ■ ■ M((3x)M(x) has strictly positive entries for some integer i > 
0. Then for any q £ M + ; the following limit exists 

cT 



1 



lim 

T^oo log T 



F(x) q dx. 



Proof. Without loss of generality, we may only consider the subsequence 
T n = (3 n . Since \G(x)\ = \F(x//3 j )\ where G is the solution of the 

associated vector equation (|1.3|), we have only to show the existence of the 
limit 

i 

lim - / \G{x)\ q dx. 

n-*co n J 

Making the change of variables x = (3 n y, we are led to prove the existence 
of the limit 

lim - ( \G{f3 n x)\ q dx. 

n— oo n J 

Notice that G(f3 n x) = M£G(x). Notice also that G(x) has strictly positive 
entries by the hypothesis on JW~. So, for the non negative matrix we 
have 

C- X \\M2\\ < \G(l3 n x)\ < C\\MZ\\ (Vx e [0, 1]) 

for some constant C > 0. By the proof of the last theorem, log \G{(5 n x)\ q dx 
is subadditive. □ 
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Remark 4.9. Let M(x) be the matrix defined by ( |1.4|) . Let M be the nu- 
merical matrix obtained by replacing fj(x) in M(x) by or 1 according to 
fj(x) = or not. Then M l > implies > 0. In particular, M d > if 
fj(x) are all strictly positive. 

Example 4.10. Let (3 > 1 be a PV-number. Let fx(x) and f%{x) be two 
strictly positive 1-periodic Holder continuous functions such that /i(0) + 
^2(0) = 1. There is a unique multiperiodic function F defined by 

™-a(sKM?)'(? 

For almost every x € M, n 1 log F((3 n x) has a limit as n — > 00; for any 
q € M + , (logT)^ 1 J Q F(x) q dx has a limit as T — > 00. 

Example 4.11. Let /5 > 1 and a, b G Z. Consider the contractive transfor- 
mations on M defined by 

_£+a _ x + 6 

D±X — - , b\X — -r= . 

For any < p < 1, there exists a unique probability measure \i with compact 
support such that 

fj, = p fx o S^ 1 + (1 — p) [i o S^ 1 . 
It is a self-similar measure. Its Fourier transform satisfies the equation 

£(s) = h(x/P)TL{x/0) + / 2 (x//3 2 )/I(x//3 2 ) 

with = £> e 27na:r and / 2 (a;) = qe 2mbx with </ = 1 — p . This is a special 

case of the equation ( |1.2| ). The corresponding matrix defined by (|1.4j ) and 
its inverse are respectively equal to 

pe 2niax qe 2nibx/f5 \ _ , / ge 2 **"/' 3 



If we take the norm |u| = max(|fi|, |f 2 1 ) on C 2 , the operator norms for M(x 



M{x) = ( f ) , M(x)" 1 = g - A e 

Dpen 

and M(x) _1 are respectively ||M(x)|| = 1 and ||M(x) _1 || = ^4?- So, when 



i- P - 

/? is a PV-number, under the condition < i, for almost all x G R the 
following limit exists and does not depend on x 

lim -log(|/I(/3 n x)| + |/I(/3 n - 1 x)|) . 

n^oo n 
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